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Abstract — The  lower  Pearl  River  in  Mississippi  and  Louisiana 
flows  through  extensive  areas  of  intertidal  marshes  and  flood- 
plains  and  is  subject  to  frequent  floods  due  to  upstream  flows, 
local  runoff  and  tidal  and  storm  forcing.  In  order  to  improve 
flow  and  flood  prediction  capabilities,  a  hydrodynamic  flow  and 
inundation  model  for  the  lower  Pearl  River  is  being  developed 
using  ADCIRC  by  the  river  modeling  group  at  the  Naval 
Research  Laboratory  at  Stennis  Space  Center.  Image  analysis 
techniques  are  applied  to  high  resolution  aerial  photography 
of  the  Pearl  River  to  extract  the  channel  boundaries  and  a 
synthetic  bathymetry  is  generated  by  fltting  an  idealized  cross- 
section  across  the  channel.  A  mesh  is  generated  for  the  entire 
lower  Pearl  River  from  Bogalusa,  LA  to  the  coast. 

A  separate  mesh  is  created  for  a  subregion  of  the  East  Pearl 
River  and  this  model  is  validated  against  observations. 

1.  Introduction 

The  Pearl  River  drains  approximately  8,600  sq.mi.  (22,270 
km^)  [  1  ]  in  southeastern  Louisiana  and  southwestern  Missis¬ 
sippi  in  the  Pine  Meadows  unit  of  the  East  Gulf  Coastal  Plain 
[2]  (see  Figure  1 ).  Originating  near  Edinburgh,  MS,  the  Pearl 
River  flows  south  before  splitting  into  two  distinct  channels 
west  of  Pieayune,  MS  to  form  the  East  and  West  branches  of 
the  Pearl,  which  drain  into  Lake  Borgne  and  the  Mississippi 
Sound  and  eventually  to  the  Gulf  of  Mexico.  The  East  Pearl 
river  serves  as  the  political  boundary  between  Louisiana  and 
Mississippi.  The  West  Pearl  is  the  major  stream  of  the  pair 
with  a  mean  flow  approximately  twenty  times  that  of  the  East 
Pearl  [3]. 

With,  generally,  70%  of  the  average  annual  rainfall  of 
68  inches  (1.72  m)  [4]  occurring  in  the  winter  and  early 
spring,  flows  are  highest  at  these  times,  with  generally  low 
streamflows  occurring  in  late  summer  and  early  fall.  The  area 
is  frequently  affected  by  tropical  weather,  including  tropical 
depressions,  storms  and  hurricanes,  which  can  bring  significant 
rain  (greater  than  30  inches  (76  cm)  [5])  as  well  as  storm  surge 
as  high  as  36  ft  (1 1  m)  [6].  With  a  tidal  range  of  as  much  as  1 
m  on  spring  tide  and  0.3  m  on  neap  at  the  ocean  mouth,  flows 
in  both  the  West,  and  especially  the  East,  Pearl  River  show 
a  strong  tidal  signal  with  tidal  fluxes  generally  much  greater 
than  river  outflow.  Both  channels  regularly  flood  and  drain  the 
bordering  intertidal  marshes. 

Flood  events  are  quite  common  in  both  branches  of  the 
lower  Pearl  River  and  can  be  caused  by  high  upstream 


Fig.  1.  The  main  channel  of  the  Pearl  River  showing  the  drainage  area  (map 
courtesy  of  Karl  Musser). 


discharge,  local  rainwater  runoff,  extreme  high  tides  due  to 
the  coincidence  of  high  spring  tides  with  Ekman  transport 
due  to  offshore  forcing  or  surge  tied  to  local  storm  events. 
Both  channels  of  the  lower  Pearl  are  bordered  by  extensive 
floodplains  which  can  be  significantly  inundated  during  high 
flood  events.  At  their  worst  these  flood  events  can  cause 
significant  physical  and  economic  damage  to  communities  on 
the  river  [7]. 

Currently,  river  stage  and  flooding  predictions  are  issued 
by  NOAA’s  Lower  Mississippi  River  Forecast  Center  using 
a  ID  hydraulic  routing  model  which  accounts  for  upstream 
flow,  direct  runoff  and  oceanic  forcing.  In  order  to  aid  in,  and 
perhaps  improve,  these  forecast  capabilities,  the  river  modeling 
group  at  the  Naval  Research  Laboratory  has  begun  to  develop 
a  two  dimensional,  depth-integrated,  hydrodynamic  model  of 
the  lower  Pearl  River  and  associated  floodplains  with  the  goal 
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of  accurately  modeling  in-channel  flows  as  well  as  flooding  in 
response  to  upstream,  oceanic  and  local  forcing.  This  model  is 
to  be  implemented  using  ADCIRC  with  an  unstructured  mesh 
derived  from  aerial  photography  of  the  Pearl  River  as  well  as 
measured  topography  and  synthetic  river  bathymetry. 

In  this  paper,  Section  II  discusses  work  done  in  generating 
a  mesh  of  the  channel  of  the  Pearl  River,  Section  HI  discusses 
work  done  to  validate  a  sub-region  of  this  mesh  and  Section 
IV  discusses  the  path  forward  and  planned  work. 

11.  Mesh  Generation 

The  section  of  the  lower  Pearl  River  to  be  meshed  runs 
from  just  below  Bogalusa,  LA  to  Lake  Borgne  and  the  Gulf 
of  Mexico  and  is  bounded  by  a  box  defined  by  30®9’N  to 
30°42’N  and  89^^30’ W  to  89°50’W.  In  the  absence  of  field 
measurements  to  establish  the  location  of  the  river  banks,  the 
use  of  image  analysis  techniques  applied  to  aerial,  or  satellite, 
photographs  of  a  river  to  extract  this  information  is  becoming 
common  [8]  [9].  The  area  is  fully  covered  by  Im  resolution 
infrared  Digital  Orthoimagery  Quarter  Quadrangle  (DOQQ) 
imagery  which  conforms  to  USGS  standards.  This  imagery  is 
housed  at  the  LSU  CADGIS  laboratory  and  is  freely  available 
from  the  Louisiana  Statewide  GIS  Atlas  website  (http:// 
atlas.lsu.edu/)  as  MRSID  compressed  georeferenced 
TIFF  imagery  on  the  UTM  15  NAD  83  datum. 

The  main  channel  and  major  side  channels  of  the  Pearl 
River  are  covered  by  twenty  five  of  these  quarter  quadrangle 
images,  a  typical  one  of  which  is  shown  in  Figure  2.  These 
images  were  decompressed  and  assembled  into  a  georefer¬ 
enced  mosaic  image  using  the  open  source  openEV  front  end 
to  the,  also  open  source,  Geospatial  Data  Abstraction  Library 
(GDAL).  This  allowed  them  to  be  assembled  into  a  seemless 
georeferenced  tiff,  32,838  pixels  by  76,674  pixels  at  1  m 
per  pixel,  which  occupied  7GB  of  disk  space  uncompressed. 
The  image  was  then  reduced  to  4  m  resolution,  requiring  900 
MB  of  space,  so  that  it  could  be  processed  using  a  desktop 
workstation  running  64  bit  MATLAB. 

A.  Image  Analysis 

As  the  images  were  taken  over  a  period  of  several  days,  with 
differing  lighting  and  exposure,  the  mosaic  was  first  equalized 
using  a  moving  patch  equalization  filter  which  significantly 
reduced  color  and  exposure  differences  across  the  image  and 
which  nearly  eliminated  lines  between  mosaic  tiles. 

The  imagery  is  saved  in  24  bit  RGB  (Red-Green-Blue)  color 
space  where  each  pixel  is  represented  by  a  value  for  red, 
green  and  blue  between  0  and  255.  However  since  HSI  (Hue- 
Saturation-Intenstity)  color  space  more  closely  matches  the 
way  that  human  eyes  perceive  images  [10],  it  is  advantageous 
to  convert  the  image  into  that  colorspace.  In  HSI  space  the 
vaue  of  hue,  ranging  from  0  to  255,  specifies  the  location  of 
the  pixel’s  color  on  a  defined  colorwheel  (with  blue  in  the 
range  of  160),  saturation  specifies  the  amount  of  grey  in  the 
color,  with  0  corresponding  to  pure  grey  and  255  to  pure  color, 
and  intensity  specifies  lightness,  or  luminance,  with  0  being 
black  and  255  being  bright  light. 


Fig.  2.  A  lypical  DOQQ  image  lile  of  the  Pearl  River  (c3008934_nes-20.sid). 
Image  center  is  al  30.4725°  N,  89.7794°  W 

Ten  subregions  were  isolated  in  the  image,  where  the  area 
of  interest  contained  similar  numbers  of  land  and  river  pixels. 
Following  a  procedure  similar  to  that  used  by  Tanaka  for  the 
Nanakita  River  in  Japan  [11],  histograms  were  generated  for 
each  channel  of  the  image  in  each  subregion.  These  histograms 
were  then  averaged  across  all  of  the  areas  for  each  channel  of 
the  image. 

A  clear  separation  between  the  hue  values  for  land  pixels 
and  water  pixels  can  be  seen  in  Figure  3,  with  water  pixels 
having  a  greenish-blue  hue  centered  around  a  value  of  135 
and  land  pixels  showing  more  spread  but  having  a  generally 
reddish-brown  hue  centered  around  a  value  of  60.  Assuming  a 
normal  distribution  for  the  hue  value  of  water  and  selecting  a 
criterion  such  that  water  pixels  should  have  a  hue  value  within 
two  standard  deviations  of  the  peak  water  hue  value,  we  can 
define  an  easily  applied  criteria  for  unsupervised  classification 
of  land  and  water  pixels. 

Applying  this  criteria  to  the  region  in  Figure  2  selects  all 
pixels  which  can  be  classified  as  water  in  this  image.  These 
results  are  then  passed  through  a  despeckling  filter,  which 
removes  isolated  small  pixel  groups  while  leaving  the  larger 
regions  which  are  indicative  of  flowing  water  or  large  ponds. 
Figure  4  shows  this  overlain  on  the  original  image  with  water 
pixels  indicated  in  red.  As  unsupervised  classification  is  not 
capable  of  discriminating  between  the  river  channel  and  larger 
ponds,  unconnected  areas  of  water,  such  as  the  pond  in  the 
lower  middle  of  Figure  4,  must  be  manually  removed.  Further 
manual  editing  is  required  to  connect  water  areas  masked  by 


Fig.  3.  Averaged  histogram  of  hue  values  in  the  vieinity  of  the  river  ehannel. 


Fig.  4.  Figure  2  with  water  pixels,  as  determined  by  applying  the  hue 
eriterion,  highlighted  in  red. 


bridges,  locks  and  trees  or  structures  along  the  channel  edge. 

Once  the  image  has  been  manually  cleaned  to  remove  these 
areas,  the  background  is  subtracted  and  a  Laplace  edge  finding 
filter  is  used  to  find  the  single  pixel  width  lines  which  bound 
the  water  area  of  interest.  This  is  illustrated  in  Figure  5. 


Fig.  5.  Water  edges  as  determined  by  a  Laplaee  edge  finding  filter. 


No  measured  bathymetry  is  available  for  the  majority  of  the 
domain,  except  for  a  portion  of  the  lower  East  Pearl  River  - 
as  will  be  discussed  later.  Lacking  any  measured  bathymetry, 
a  synthetic  bathymetry  must  be  generated  for  the  mesh.  For 
each  point  identified  as  being  a  water  point,  approximately  4 
million  in  the  full  river,  a  parabola  [13]  is  fit  spanning  the 
channel  along  the  shortest  possible  path  and  an  initial  depth 
is  taken  from  that  parabola.  After  the  initial  bathymetry  is 
calculated,  each  point  is  averaged  with  its  nearest  neighbors  to 
insure  smooth  transitions  and  avoid  discontinuities,  especially 
in  areas  where  channels  split  or  merge. 

This  ordered  coastline  and  synthetic  bathymetry  are  then 
used  to  generate  an  initial  mesh.  At  the  time  of  publication, 
this  initial  mesh  is  still  undergoing  mesh  refinement  and  has 
yet  to  be  run.  It  is  estimated  that  the  final  mesh  will  have 
approximately  600,000  nodes  and  a  mesh  resolution  which 
ranges  from  4  m  in  the  smallest  channels  and  side  creeks  to 
500  m  at  the  ocean  boundary. 


B.  Mesh  Generation 

In  order  to  generate  the  unstructured  mesh  we  employ  the 
NRL  Mesh  Create  program  [12]  which  requires  an  input  of 
an  ordered  coastline  as  well  as  bathymetry  which  spans  the 
entire  domain.  An  ordered  coastline  can  be  extracted  from  the 
raster  edge  data,  such  as  is  seen  in  Figure  5,  by  employing  an 
edge  following  algorithm  to  trace  each  edge  and  output  it  as 
a  clockwise  or  counterclockwise  ordered  line. 


111.  Validation 

As  the  full  mesh  is  still  being  developed,  as  of  the  time  of 
this  writing,  a  separate  sub-mesh  was  selected  for  validation 
of  the  model  and  method.  This  mesh  covers  the  portion  of  the 
East  Pearl  River  near  and  downstream  of  Stennis  Space  Center 
and  is  in  a  region  for  which  detailed  bathymetry  is  available 
along  with  some  amount  of  validation  data. 


Topography  (♦)  /  Bathymetry  (•) 


£  30  28 

I  30  26  .5 

g  30  24 
30  22 
302 

30  Id  H 

3016  H 

^  ■  15 

■09  7  89  65  -89  6  89  55  89  5 

l  ongrtixJo  ( W) 

Fig.  6.  Measured  bathymetry  in  the  East  Pearl  River  subregion. 


A.  East  Peart  River  Sub-mesh 

Detailed  bathymetry  is  available,  from  a  number  of  sources, 
for  the  lower  portion  of  the  East  Pearl  River  from  just  upstream 
of  Stennis  Space  Center  to  the  ocean  boundary  (see  Figure  6). 
Where  bathymetric  data  sets  overlap  they  have  been  averaged 
to  create  a  single  bathymetry  file  covering  the  whole  of  the 
lower  East  Pearl  region.  This  has  been  combined  with  edge 
data  from  a  mesh  of  the  area  previously  created  at  NRL 
to  generate  a  new  mesh  of  the  area  of  interest.  This  is  an 
unstructured  mesh  with  80,452  nodes  and  143,645  elements. 
Grid  resolution  varies  from  approximately  1  m  in  the  smaller 
channels  to  500  m  in  Lake  Borgne  and  in  the  Mississippi 
Sound. 

The  upstream,  or  river  flux,  boundary  is  modeled  using  a 
defined  flux  (IBTYPE  =  2)  boundary  condition  which  imposes 
a  defined  normal  flow  at  the  boundary  nodes.  As  this  has  the 
effect  of  strongly  damping  the  along  channel  tidal  flows,  which 
are  not  able  to  radiate  through,  and  causes  the  flooding  tidal 
waver  to  reflect  and  interfere  with  predicted  elevations  through 
the  domain,  it  is  necessary  to  artificially  move  this  boundary 
far  from  the  region  of  interest  such  that  this  damping  effect 
becomes  minimal  in  the  true  model  domain  [14].  For  this 
reason  the  mesh  was  extended  upstream  as  a  virtual  channel 
for  a  distance  equal  to  1.5  times  the  main  channel  length.  This 
distance  was  found  to  be  sufficient  such  that  tidal  flows  are 
unimpeded  through  the  river  domain  and  the  reflected  tidal 
wave  is  dissipated  before  it  interferes  with  elevations  in  the 
region  of  interest.  The  mesh,  with  this  extension,  is  shown  in 
Figure  7.  River  forcing  is  set  to  a  constant  value  of  5(X)  efs. 
This  corresponds  to  the  mean  measured  flow  near  the  upstream 
boundary  of  the  mesh  for  the  last  two  weeks  of  May  2(X)9. 

Tidal  forcing  is  applied  to  the  downstream  open  boundaries 
with  six  component  harmonic  tides  (Ql,  01,  Kl,  N2,  M2, 
K2)  taken  from  the  Grenoble  FES99  ADCIRC  database.  Many 
of  the  small  side  channels,  particularly  in  the  lower  portion 
of  the  mesh,  feed  into  small  ponds  and  marshes.  As  such 
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Fig.  7.  Easl  Pearl  River  submesh  wilh  upslream  exten.sion. 

these  have  been  left  open  as  radiation  boundary  conditions 
allowing  the  tides  to  radiate  through  them.  The  period  of 
interest,  18-31  May  2009,  was  selected  as  a  period  during 
which  meteorological  forcing  was  minimal  and,  as  such,  it 
has  not  been  included  in  this  test  run. 

B.  Data  Sources 

Data  to  validate  the  model  is  available  from  a  number 
of  sources  spanning  the  length  of  the  main  channel  (see 
Figure  8).  NRL  maintains  a  long  term  flow  mooring  near  the 
upper  boundary  of  the  mesh  consisting  of  a  side  looking  RDl 
HADCP  which  provides  instantaneous  cross-channel  averaged 
velocities  and  estimated  volume  flux.  The  USGS  maintains  a 
streamflow  monitoring  station  (number  301 14 10893 203(X))  at 
the  CSX  railroad  bridge  near  Clairbome,  MS  which  provides 
water  depth  and  streamflow  near  the  ocean  boundary.  During 
the  months  of  May  and  June  2(X)9  two  additional  moorings 
were  placed  along  the  axis  of  the  East  Pearl  River  (at  the  1-10 
and  Highway  90  bridges).  These  bottom  mounted  Aquatroll 
CTDs  recorded  water  depth  every  30  minutes  throughout 
their  deployment.  There  is,  additionally,  a  USGS  streamflow 
monitoring  station  located  near  the  HADCP  mooring  which 
intermittently  provides  water  depth  and  streamflow  informa¬ 
tion. 

C.  Validation 

The  model  was  run  in  2D  depth-integrated  barotropie  mode 
for  the  period  from  18  to  31  May  2009.  River  flux  and  tides 
were  ramped  up  over  a  21  day  period  prior  to  the  period 
of  interest.  Velocity  and  depth  were  output  for  each  node 
every  hour.  As  the  model  run  did  not  account  for  subtidal 
variations  tied  to  forcing  by  wind  and  atmospheric  pressure. 


Fig.  8.  Locations  of  East  Pearl  River  moorings  during  May  and  June  2009. 
Figure  courtesy  of  Wesley  Goode. 


all  measured  data  was  filtered,  using  a  3rd  order  40  hour  high 
pass  Butterworth  filter,  to  remove  subtidal  signals  [15].  Once 
these  subtidal  signals  were  removed,  the  model  and  the  data 
show  very  good  agreement  for  both  tidal  variations  in  water 
depth  and  along  ehannel  veloeity. 

Figure  9  shows  tidal  variation  in  water  depth,  referenced 
to  mean  water  level,  at  the  CSX  railroad  bridge  from  18-31 
May  2009.  The  blue  curve  is  the  measured  data  and  the  red 
curve  the  model  output.  The  model  matches  the  phase  of  the 
measurements  preeisely  while  overpredieting  the  magnitude 
of  the  spring  tide  by  approximately  3  %  and  matching  it  very 
well  at  neap  tide. 

Further  upstream,  at  the  1-10  Aquatroll  mooring,  we  see 
similar  results.  Figure  10  shows  the  same  data  plotted  for  the 
1-10  mooring.  As  at  the  lower  boundary,  the  model  slightly 
overprediets  tidal  height  at  spring  tide,  in  this  ease  by  approx¬ 
imately  5  %.  At  neap  tide  the  high  tide  is  modeled  preeisely 
but  the  low  tide  is  overpredieted,  again  by  approximately  5  %. 
Both  other  mooring  loeations  show  similar  results. 

Comparing  the  measured  eross  ehannel  average  of  the  along 
ehannel  veloeity  at  the  HADCP  mooring  to  the  model  output 
(see  Figure  1 1 )  the  agreement  is  also  good,  however  the  model 
underprediets  maximum  ebb  (positive)  veloeities  by  as  much 
as  15  %  on  spring  tides.  This  is  possibly  due  to  inaeeuraeies 
in  modeling  river  flux  as  eonstant  in  this  case  or  perhaps 
to  surface  effects  interfering  with  the  eross  ehannel  HADCP 
veloeity  measurements. 

Current  maps  produced  for  different  regions  of  the  Pearl 
River  (see,  for  example,  Figures  12  and  13,  the  locations  of 
which  are  shown  in  Figure  14)  show  the  expected  flow  patterns 
but  there  is,  at  present,  no  data  available  to  validate  these. 

D.  Validation  of  the  Synthetic  Bathymetry 

In  order  to  test  the  validity  of  the  simple  synthetic 
bathymetry  which  has  been  used  in  the  main  mesh,  the 
real  bathymetry  of  the  validation  model  was  replaced,  node 


Fig.  9.  Tidal  variations  in  waier  depth  at  the  CSX  railroad  bridge  from  18 
10  31  May  2009. 
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Fig.  10.  Tidal  variations  in  water  deplh  at  ihe  I- 10  mooring  from  18  to  31 
May  2009. 


0.12 


Fig.  11.  Cross  channel  average  of  the  along  channel  velocity  at  the  HADCP  pig.  13.  Current  naap  in  the  middle  region  of  the  East  Pearl  mesh  at  1200 


mooring  from  18  to  31  May  2009. 


CST,  20  May  2009. 


J.:  iOn-^  '•■in'' 


Fig.  12.  Current  map  in  the  upper  region  of  the  East  Pearl  mesh  at  1200 
CST.  20  May  2009. 


Fig.  14.  Location  of  current  maps  in  Figures  12  and  13. 
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Fig.  15.  Topographically  adjusted  synthetic  bathymetry  for  the  East  Pearl 
River  mesh. 

by  node,  with  a  bathymetry  based  on  a  parabolic  channel 
cross-section  as  described  in  Section  11-B.  As  this  synthetic 
bathymetry  is  referenced  to  mean  sea  level  it  docs  not  account 
for  the  topographic  slope  of  the  area  and  thus  does  not  allow 
for  modeling  of  the  gravity  driven  pressure  gradient  attributed 
to  that  slope.  To  correct  for  this,  the  synthetic  bathymetry 
was  merged  with  topographic  data  from  the  NGLl  (Northern 
Gulf  Littoral  Initiative)  three  are  second  topography  of  the 
region.  The  edges  of  the  parabolic  cross-sections  were  fixed 
to  the  topographic  height  from  the  NGLI  data  and  the  channel 
was  modeled  downwards  from  this  new  datum.  This  synthetic 
bathymetry  is  shown  in  Figure  15. 

Comparing  this  synthetic  bathymetry  with  the  measured 
bathymetry  shown  in  Figure  6,  it  appears  that  the  syn¬ 
thetic  bathymetry  is  a  good  approximation  of  main  channel 
bathymetry,  ineluding  the  scour  holes  near  channel  junctions, 
but  that  it  does  not,  and  in  fact  ean  not,  model  either  the 
shoaling  regions  in  the  channel  or  the  ehannel  meanders  in 
whieh  the  thalweg  is  displaced  to  the  left  or  right  of  channel 
center.  The  thalweg  is  constrained  to  be  in  the  geometric  center 
of  the  channel. 

At  the  ocean  boundary,  the  synthetie  bathymetry  does  not 
model  the  deep  scour  channel  created  by  the  high  tidal  flows 
in  the  area  and  it  significantly  overestimates  water  depth, 
especially  in  Lake  Borgne.  As  this  bathymetry  is  based  purely 
on  horizontal  geometry  it  can  not  model  the  sedimentation  and 
scour  processes  whieh  reduce  or  increase  water  depth  in  these 
areas.  However  since  these  are  in  regions  of  lesser  interest  in 
this  modeling  effort,  and  are  near  boundaries  which  are  forced 
by  tidal  elevation  referenced  to  mean  sea  level,  there  should 
be  little  effect  on  the  overall  solution. 


It  is  apparent,  especially  near  the  upstream  boundary,  that 
the  synthetic  bathymetry  overestimates  the  topographic  height. 
This  is  likely  due  to  the  low  resolution  of  the  three  arc  second 
NGLl  topography.  However  this  error  amounts  to  only  2.5  m 
over  a  length  from  the  coast  of  24.5  km  and  is  not  likely  to 
be  significant.  For  future  work  we  are  examining  both  USGS 
topography  and  LIDAR  derived  digital  elevation  models  of  the 
area  as  well. 

Removing  the  topography  from  both  sets  of  bathymetry, 
such  that  both  are  referenced  to  mean  sea  level  at  all  cross- 
scctions,  and  examining  only  the  region  of  the  main  channel 
upstream  of  the  obvious  differences  due  to  tidal  and  sediment 
processes,  there  is  a  maximum  error  in  the  main  channel  depth 
of  0.72  m  where  the  synthetic  bathymetry  does  not  account 
for  a  shoal  in  the  river.  The  RMS  error  for  all  points  in  the 
main  channel  is  0.21  m. 

Near  the  ocean  boundary  and  in  Lake  Borgne  there  is  a 
maximum  error  of  10.13  m  where  the  synthetic  bathymetry 
greatly  overpredicts  depth  in  the  shallow  and  sediment  filled 
Lake  Borgne. 

IV.  Future  Work 

Final  refinement  of  the  mesh,  including  both  the  East  and 
West  branches  of  the  lower  Pearl  River,  is  expected  to  be 
complete  shortly.  The  synthetic  bathymetry  generated  for  this 
mesh  will  be  offset  by  the  mean  topography,  as  was  described 
in  Section  Ill-D,  and  then  merged  with  the  available  measured 
bathymetry  for  the  lower  seetion  of  the  East  Pearl  River,  shown 
in  Figure  6.. 

This  new  model  will  then  be  validated  against  the  mea¬ 
surements  taken  by  the  permanent  moorings  in  the  East  Pearl 
River  (described  in  Section  III-B)  and  planned  measurements 
of  water  depth  to  be  taken  in  the  heretofore  ungauged  West 
Pearl  River. 

Following  this  validation,  the  mesh  will  be  extended  to  in¬ 
corporate  the  floodplain  around  and  between  the  two  branches 
of  the  river  using  both  USGS  and  NGLl  topography.  The 
2005  NOAA  landcover  database  provides  information  on  the 
dominant  local  vegetation  which  can  be  correlated  to  drag 
coefficients  [16]  [17]. 

Once  this  work  is  complete,  higher  flowrates  will  be  mod¬ 
eled  such  that  the  river  overflows  its  banks  and  modeled 
inundation  of  the  floodplain  will  be  compared  to  historical 
records  of  flooding. 

Work  is  underway  to  develop  the  model  to  run  in  real  time, 
providing  daily  updated  seven  day  forecasts  of  river  stage, 
tidal  elevation,  currents  and  inundation  which  will  be  made 
publicly  accessible  on  the  web.  This  information  can  then  be 
used  to  better  constrain  and  provide  boundary  forcing  to  the 
NOAA  model  currently  in  use  to  forecast  for  the  West  Pearl 
River. 

V.  Conclusion 

A  methodology  has  been  developed  allowing  a  mesh  of 
a  large  region  of  the  lower  Pearl  River  in  Mississippi  and 
Louisiana  to  be  developed  at  4  m  resolution  from  aerial 


photographs.  In  the  absence  of  measured  bathymetry  a  “best 
guess”  synthetic  bathymetry  has  been  fit  to  the  channel  by 
assuming  a  standard  parabolic  profile. 

A  subregion  of  this  mesh  has  been  validated  against  existing 
measurements  and  been  shown  to  accurately  predict  water 
depth  and  currents  at  all  measurement  locations.  The  method¬ 
ology  to  generate  a  synthetic  bathymetry  has  been  validated 
against  the  measured  bathymetry  in  the  East  Pearl  River  and 
shown  to  be  a  realistic  model  in  the  main  river  channel  but  to 
be  very  inaccurate  in  coastal  regions  where  sedimentation  and 
scour  can  significantly  affect  channel  morphology. 

Plans  arc  in  place  to  further  develop  this  mesh  and  extend 
it  into  the  neighboring  floodplains  to  allow  modeling  of  the 
inundation  of  these  plains  tied  to  river  flux,  local  runoff  and 
storm  events. 
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